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O ! Abstract 

^! Given n polynomials in n variables with a finite number of complex roots, for any of their 
,Jh^ roots there is a local residue operator assigning a complex number to any polynomial. This 
is an algebraic, but generally not rational, function of the coefficients. On the other hand, 



> 



the global residue, which is defined as the sum of the local residues over all roots, has 
invariance properties which guarantee its rational dependence on the coefficients [9], [27]. 
In this paper we present symbolic algorithms for evaluating that rational function. 

Under the assumption that the deformation to the initial forms is flat, for some choice of 
weights on the variables, we express the global residue as a single residue integral with 
respect to the initial forms. When the input equations are a Grobner basis with respect to 
a term order, this leads to an efficient series expansion algorithm for global residues, and 
to a vanishing theorem with respect to the corresponding cone in the Grobner fan. 

The global residue of a polynomial is shown to equal the highest coefficient of its (Grobner 
basis) normal form, and, conversely, the entire normal form is expressed in terms of global 
residues. This yields a new method for evaluating traces over zero-dimensional complete 
intersections. Applications to be discussed include the counting of real roots (as in [4], [22]), 
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the computation of the degree of a polynomial map (cf. [12]), and the evaluation of mul- 
tivariate symmetric functions (cf. [16], [21]). All results and algorithms are illustrated for 
an explicit system in three variables. 
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0. Basic Properties of Multidimensional Residues 



Multidimensional residues play a fundamental role in complex analysis and geometry. Re- 
cent applications of residues in computational algebra include explicit division formulae 
[5], [6], the evaluation of symmetric functions [21], the membership problem for polyno- 
mial ideals [9], [10], the effective Nullstellensatz [13], and numerical algorithms for solving 
polynomial systems [7]. In most of these articles the emphasis lies on degree estimates 
and complexity results. Our goal here is to develop practical tools for computing global 
residues. We thus refrain from using "univariate projections" or "linear changes of co- 
ordinates"; instead we seek algorithms involving Grobner bases and sparsity-preserving 
series expansions. While initially our discussion follows a path similar to [25], [27], it then 
proceeds to systematically develop the interplay between residues and Grobner bases. 

This paper is organized as follows. In §1 we consider n polynomials in n variables which 
form a Grobner basis (or H-basis) with respect to some choice of positive weights. This 
hypothesis has natural geometric (1.3'), algebraic (1.5) and analytic (1.7) interpretations. 
In (1.17) we express the global residue as a single residue integral with respect to the initial 
forms. In §2 we specialize to the case of a Grobner basis in the usual sense, with respect 
to a term order. In (2.3) we express the residue as a coefficient of a certain polynomial. 
This yields a polyhedral vanishing theorem (2.5), and a bound on the degree of the residue 
as a polynomial in the trailing coefficients (2.7). The constructions of §1 and §2 lead to 
algorithms, which will be presented in §3. In §4 we relate global residues to the coefficients 
in the (Grobner basis) normal form. The global residue of a polynomial is shown to 
equal the highest coefficient of its normal form (4.2). This results in fast procedures for 
computing residues and traces (4.8). In §5 we present applications to symmetric functions, 
to computing the degree of a polynomial map, and to counting real roots. Finally, we 
study in §6 an explicit system in three variables. 

In this section (§0) we review the complex-analytic definition and basic properties of mul- 
tidimensional residues. Details and proofs can be found in [1],[15],[25]. For the algebraic 
counterpart to the analytic theory see e.g. [3], [19], [20], [23]. We shall discuss the equivalence 
of the algebraic and the analytic approach briefly at the end of §0. 

Given n holomorphic functions gi,...,g n in an open set U C C n with a single common 
zero p in U, one can associate to any holomorphic function h G 0(U) the local residue at 
p of the meromorphic n-form 



uj = 



h(x) dx. 



0l(x)...0 n (x) ' 



<ix = dx\ A ... A dx r 



This defines the C-linear operator 

(0.1) 0{U)^C, h h-> Res p ( h W dx \ := _L_ f u 

where r g (e) is the real n-dimensional cycle 

r g( e ) = {xGt/" : |&(x)| = e< ,i = l,...,n} 

with orientation defined by the n-form da,rg(gi) A. . . Aciarg(gr n ), and where e = (ei, . . . , e n ) 
is any n-tuple of sufficiently small, positive real numbers. 

The Cauchy formula in n complex variables provides the simplest example of a residue 
operator: If gi(x) = (xi — Pi) ai+1 , ai G N, then 

, nn . „ / h(x)dx \ 1 /^i+.-.+^n 

(a2) ReSp J = Ur-feg- J (p) ■ 

Let J g := det denote the Jacobian of g = (<7i, . . . ,g n )- Then, 

(0.3) Res p ( feJgdx ) = // g (p) %) , 

V^l ■■■9nJ 

where jtt g (p) denotes the intersection multiplicity of g = (g±, . . . , g n ) at p [15, p. 662ff.]. If 
p is a simple root, hence J s (p) ^ 0, then 

/ \ _ / hdx \ hip) 

(0.4) Res p ' 1 - ' 



gi---g n J J s (p)' 

Suppose now that gi,...,g n G C[x] are n-variate polynomials whose zero set Z(g) is a 
non-empty finite subset of C n . We can consider the global residue operator ([9], [15], [25]) 
which assigns to a polynomial h G C[x] the complex number 

^ / . \ „ f hdx \ \ ~\ „ f hdx \ 
(0.5) Res,,,) = Res (— ) == Re Sp (_) . 

The main functorial properties of the global residue are encapsulated in the following two 
results whose proof may be found, for example, in [25, II. 8. 3-4]. 

(0.6) Transformation Law. Suppose that fi,...,f n G C[x] have finitely many common 
roots and that we can write 

n 
3 = 1 

Then, for h G C[x], 

„ f hdx \ _ f h det(An) dx 
Res = Res ' V j! 



g\---g n ) \ h---f; 
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If <7i, . . . ,g n G C[x] are as above, the ideal 7 generated by them is zero-dimensional, and 
therefore V = C[x]/7 is a finite dimensional C- vector space. Since the global residue 
Res g (/i) vanishes for h E I (see [15, p. 650]), it defines a C-linear map: 



Res s :V -> C , Res s (h) . 



(0.7) Duality. A polynomial h G C[x] lies in 7 if and oniy ifRes s (fh) = for ail / G C[x] . 
This duality law may be interpreted as follows: Let V* = Horned C). Then the pairing 



where (6.0) (6') = 0(66'), makes V* into a F-module. Statement (0.7) is equivalent to the 
assertion that the residue operator Res g is a generator of V* as a F-module. 

We recall that Res g (6) is a rational function, with integral coefficients, in the coefficients 
of gi,...,g n . It may be computed in simply exponential time with respect to n, the 
number of variables; see [9], [21], [27]. A general procedure, suggested in these articles, 
is to find univariate polynomials fi(xi), . . . , f n (x n ) in the ideal generated by gi, ■ ■ ■ ,g n 
and to transform the global residue applying (0.6). The global residues with respect to 
/i,...,/n are then computed as a sum of products of univariate global residues. This 
general procedure is often too slow for practical computations. We seek more efficient 
algorithms for "nice" situations, such as the case when gi,...,g n are a Grobner basis for 
a term order. The theory for such nice situations is to be developed in the next section. 

In closing let us mention the relationship between this analytic definition of global residue 
and the algebraic definitions: of course, they coincide. Write for each i = 1, . . . , n, 



and denote A := det(^y). Let U C C n be the union of relatively compact open neigh- 
borhoods isolating each of the points in Z(g), let e be any n-tuple of small positive real 
numbers, and define II e := {x G U : |<7i(x)| < e\ ,Vz = 1, ...,n}. For any holomorphic 
function h on U, one can deduce from (0.2) and (0.6) (cf.[25, §17]) the following integral 
representation known as Weil's formula [26]: 



VxV* 



V 



* 



(6, 0) - (6.0) , 



n 
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As |<7i(x)| < |(7i(y)| for any i , x 6 II £ and y G r g (e), the integrand may be expanded as a 
multiple geometric series 

/i(y)A(y,x) \X^fTT #*( x ) Qi 

n^w-aw) = ' l(y)A(y ' x) S(S^ TT 

which converges uniformly on compact subsets of r g (e) x n e . As a result of term-by-term 
integration, we deduce that (cf.[25,§5],[5],[6]) 



(0.8) fc(x) = Res s (h( ■ ) A( • ,x)) + £ Res f fft^-^ ) g C 

iai>i viii=ifAyj / 



(x), vxen £ . 



Note that the second summand on the right is in the ideal generated by g\, . . . ,g n in 0(H € ) 
and that Res g (/i( • ) A( ■ , x)) depends polynomially on x. If, in addition, h is a polynomial, 
then the fact that Z(g) is contained in II e plus the fact that local analytic membership is 
equivalent to local algebraic membership ([24]), imply that 

(0.9) /i(x) = Res g (/i( • ) A( ■, x)) on the quotient ring V . 

In general, (0.8) does not provide a representation of their difference as a polynomial linear 
combination of gi, . . . , g n . Under the hypothesis (1.3) below, it follows from the vanishing 
statement in (1.18) that the series becomes a finite sum, giving an effective division formula 
with remainder which involves computing only finitely many global residues associated to 
powers of gi, . . . ,g n . In summary, the formula (0.9) is the algebraic version of the integral 
representation. It proves that the global residue we are considering coincides with the 
"trace" associated to A as in [19, Appendix F] and [13], and with the Kronecker symbol 
(i.e., the dualizing linear form associated to A) as in [3]. 



1. Grobner Bases for a Weight Partial Order 

Let K be any subfield of the complex numbers C and let <7i(x), i = 1, . . . , n, be polynomials 
in S = K[x], x = (xi, . . . , x n ) . Let w = (w\, . . . , w n ) G N n be a positive weight vector. 
The weighted degree of a monomial x a = x^ 1 . . . x® n is 

n 

deg w (x a ) = (w,a) = ^Wim. 

i=i 

We extend the notion of weighted degree to arbitrary polynomials in S in the usual manner. 
Write each polynomial #i(x) as 

(1-1) £i(x) =Pi(x) + 9i(x) 
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where pi is w-homogeneous, and 



(1-2) di = deg w (p l ) = deg w (^) ; deg w (^) < deg w (^) . 

Throughout this section we make the following assumption: 

(1.3) = • • • = p n (x) =0 if and only if x = 

In what follows we will interpret this condition geometrically (1.3'), algebraically (1.5) and 
analytically (1.7). Let 

(1.4) ut-^)=t dl g l {t~ wi xi,...,t- w -x n ) 

This is a homogeneous polynomial of degree di in (t;xi,...,x n ) relative to the weights 
(l;w!, . . . ,w n ). 

Let denote the weighted projective space with homogeneous coordinates (t; 

and weights (l;wi, . . . , w n ). The image of the hyperplane {t = 1} C C n+1 \{0} in P™ is 

identified with C n . If 

A = {(t;x)eP; : &(*;x) = 0} 
then (1.3) is equivalent to the geometric condition 

(1.3') Din...n^cC 

The algebraic meaning of (1.3) is best expressed using the following notion of a Grobner 
basis: Given a polynomial / e S, we denote by in w (/) its form of highest weighted degree. 
For any ideal I C S we define the initial ideal in w (I) to be the ideal generated by in w (/) 
where / runs over /. A finite subset Q C / is said to be a Grobner basis for /, relative to 
the weight w, provided: 

in w (J) = (m w (g) : g e Q) . 

We emphasize that in w (J) need not be a monomial ideal. Some authors prefer to call Q 
an H-basis, a term which goes back to Macaulay in the classical case w = (1, 1, . . . , 1). 

(1.5) Lemma. Suppose Q = {gi, . . . ,g n } C S satisfy (1-3). Then Q is a Grobner basis 
for the ideal it generates. Conversely, suppose Q = {gi, . . . , g n } C S is a Grobner basis , 
with respect to w, for a zero-dimensional ideal I. Then {gi, . . . , g n } satisfy (1.3). 
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Proof: With the same notation as above we have pi = in w (^j) and qi = gi — in w (^j). Since 
Pi, • • . ,p n define a complete intersection, the Koszul complex on these forms in exact. This 
implies that every syzygy Yl^i' e i on (Pi> • • • > Pn) can De written as a linear combination 
of the basic syzygies pj • et — Pk • £j- 

Suppose that Q is not a Grobner basis. Then, there exists a polynomial 

/ = S ^ j h i g i = ^hipi + ^h i q l 

whose initial form does not lie in (pi, . . . ,p n ). Hence Y2^iPi = 0- By the remark above, 
we can write ^2 hi ■ ei = J2j k ^jk • (Pj • e& — Pk • and the leading term of 

^higi = ^hiqi = ^ b jk • (PjQk ~ PkQj) 

i,k 

must lie in (pi, . . . ,p n ). This is a contradiction, completing the proof of the first statement. 
To prove the converse, it suffices to note dim (J) = dim(in w (J)) = (see e.g. [17]). o 

(1.6) Remarks: (i) The first part of Lemma (1.5) remains true for any set of polynomials 
G = {<7i) • • • ? dm} C S whose initial forms define a complete intersection. 

(ii) The initial ideal in w (7) is a flat deformation of the given ideal I (see e.g. [11, Ch. 6]). 

(iii) The results in this section can be extended to fields other than the complex numbers 
using the deformation techniques in [20]. 

For each t e C, consider the map g t : C n — > C n defined by gt(x) = (gi(t; x), . . . , g n (t; x)). 
We have the following analytic interpretation of (1.3). Recall that a map F:C n — > C n is 
said to be proper if the inverse image of any compact set is compact. 

(1.7) Lemma. The polynomials gi, . . . , g n satisfy condition (1.3) if and only if the map 
gt is proper for every t £ C. 

Proof: At t = we have go = (pi, . . . ,p n ). Since the polynomials pi are weighted 
homogeneous, the inverse image g ~ 1 (0) is compact if and only if g ~ 1 (0) = {0}. Thus, if 
go is proper, then condition (1.3) is satisfied. 

For the converse it is enough to show that the map g is proper, since (1.3) is a condition 
on just the initial form of the polynomials. Let g: C n+1 \ {0} — > C n+1 \ {0} be defined by 

H(t,x 1 ,...,x n ) = (t,£i(t;x),...,£ n (t;x)) 
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The fact that g(x) satisfies (1.3) guarantees that g takes values in C n+1 \ {0}. Since 



g defines a map from P™ to weighted projective space with weights (1; d±, . . . , d n ). We 
may now consider the embedding of C n in C n+1 \ {0} as the hyperplane {t = 1}. Since 
t has weight one in both P™ and P^, the natural projection from C n+1 \ {0} to P™ or 
P^ is a homeomorphism of the hyperplane {t = 1} to its image. Thus, g is a continuous 
extension of g to appropriate compactifications of C n . If K C C n is compact then g _1 (if ) 
is compact since it coincides with g -1 (if). o 

(1.8) Examples: An important special case, to be investigated in detail in §2, is that of 
n polynomials gi, ■ ■ ■ ,g n £ S with finitely many common roots in C n and such that they 
are a Grobner basis with respect to some term order -<. We can choose a weight vector 
w G N n such that in^(gi) = in w ((7i). Since the ideal generated by the initial monomials 
in w ((7j) is zero-dimensional, we may assume without loss of generality that: 

m w (gi) = aixY +1 

for some ^ G K \ {0}, and therefore they satisfy (1.3). Particular examples are: 

(i) The term order -< is lexicographic order: then gi, ■ ■ ■ ,g n satisfy #i(x) = gi(xi, . . . , x n ) 
and gi is monic in Xi. This is the case studied in [9]; see also [21] for the subcase when 
<7i(x) = gi(xi) are univariate polynomials. 

(ii) The term order is defined as total degree with ties broken by lexicographic order with 
x n > . . . > x\. This is the case studied in [1, (21.3)] and [25, II. 8. 2]. A weighted variant, 
due to Aizenberg and Tsikh, is studied in [1, (21.5)]. 

In this section we are interested in studying the global residue Res g (/i), for a polynomial 
h G C[x], under the hypothesis that (?i(x), . . . , g n (x) satisfy (1.3). This hypothesis makes 
it possible to reduce the computation of Res g (/i) to that of residues involving only certain 
powers of the initial forms £>i(x), . . . ,p n (x). This is the content of (1.20) below. In fact, 
considering t as a parameter, the idea that one can recover the information from the 
deformation to the initial forms is the core of the geometric interpretation of Grobner 
bases (see e.g. [2]). 

Since the map g(x) = (gi (x), . . . , ^ n (x)) is proper, we can replace (0.5) by a single integral 



g(At, X Wl xi,..., X w - x n ) = (At, A dl hit; x), • • • , A d " g n (t; x)) , 
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for any r = (r 1? . . . , r n ) G (M >0 ) n , where T(r) is the compact, real n-cycle 

r(r) = {x£C" : |<7i(x)|=r i; i = l,...n}. 

Similarly, for each fixed t G C, 

h(x) <ix 



(1.9) R h (t):=Re S ( , , \ = * [ 



9n(t;x)J (2Tvi) n Jr t(r) sfi(t;x)...£ n (t;x) ' 

where the global residue is taken relative to the divisors { x G C n : gi(t;x) = 0}, z = 
1, . . .,n and f t (r) = { x G C n : |&(*;x)| = n ; i = 1, . . .n}. 

The family {(?i(£, x), . . . , g n (t, x)} is a Grobner basis with respect to w' = (1, 2w±, . . . , 2w n ) 
for the ideal I generated by {/, / G /} and in w /(/) = in w (/), V/ G 5. It then follows that 
the coordinates (£, x) are in Noether position [8] for 7 and, consequently, we may apply 
Theorem 3.3 in [9] to deduce that Rh(t) is a polynomial in t. We will reprove this and in 
fact obtain the stronger result (1.18). 

We set 

n d w 

(1.10) G(f;x) := n^(*! x ) = Il A j(x)t J 

i=l j=0 

where d w = d\ + . . . + <i n . Inverting the polynomial G(t; x) as a rational formal power 
series in t, we write 

(1.11) G-^fjx) = ^S,(x)^. 

J>0 

Given positive real numbers fci, . . . , fc n , let 

T(k) :={xGC" : | Pi (x) | = h } . 



(1.12) Lemma. Given 5 > 0, there exist positive constants k±, . . . , k n so that, for \t\ < 5, 
and this series is uniformly convergent for \t\ < 5. 
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Proof: Suppose gi, ■ ■ ■ ,g n satisfy (1.3) and let 

&(*;x) =Pi(x) +t^(t;x). 
Then, given 5 > 0, there exist positive constants fci, . . . , k n such that, for x e T(k), 

(1.13) ^>N l (t;x)| 

for alH = 1, . . . , n and \t\ < 5. Indeed, because of the weighted homogeneity property of 
Pi, it suffices to take ki = X di with A sufficiently large. 

The estimate (1.13) allows us to apply Rouche's principle for residues [25, II. 8.1] and 
replace, for \t\ < 5, the integration cycles r t (r) in (1.9) by the fixed cycle T(k). Thus, 

(1.14) R h (t) = — ^— [ , h ^ dx f or a n \ t \ < s. 

1 ' K> (2ttz)- 7 T(k) gi (t-x)...g n (t;x) 1 1 " 

In view of (1.13), it follows that the series 

V^-(x)^ = j± : 

is uniformly convergent for x G T(k) and \t\ < 5. Since we can now integrate (1.14) 
term-by-term, the result follows, o 

(1.15) Lemma. Let -P(x) = pi(x) . . .p n (x). Then P m+1 (x) S m (x) is a weighted homo- 
geneous polynomial of degree m (<i w — 1) with respect to w. 

Proof: Since G(t; x) is weighted homogeneous of degree <i w and t has weight 1, the 
coefficients A,(x) in (1.10) are weighted homogeneous of degree d w — j. On the other 
hand, the series (1.11) inverts (1.10). This implies the following recursion relations: 

m 

(1.16) ^2 A J B m-j =0 , TO > 1 

i=o 

with initial conditions Aq Bq = 1 and Ao(x) = -P(x). In particular, 

m 

and 

P™ +1 B m = ^A J P^- 1 (P^+ 1 J B m _ J ). 
j'=i 

Assuming that (1.14) holds inductively with respect to to, we obtain 

deg w (P m+1 S m ) = (rf w -j) + (j-l)c/ w + (m-j)(rfw-l) = m(dw-l). o 

The following is the main result in this section. 
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(1.17) Theorem. For any monomial x a = x± x . . . x°^ , set s(a) = (w, a) — d w + Ym=i w i- 
Then 

R^(t) = Res g (x a ) -t s ^ , 



and 

Res g (x a ) = — — — / x a B s(a) (x)dx. 
{2m) n y T(k) 



Before proving (1.17), we note the following weighted version of the Euler-Jacobi Theorem 
[15, p. 671]. A more general toric version was given by Khovanskii in [18]. 

(1.18) Corollary. Rh(t) is a polynomial in t of degree at most deg w (/i) — d w + Y Wi, and 

Res ( — — J = whenever deg w (/i) < d w — Y\ w i- 

\9i ■ ■■9nJ 

We observe also that, under the current hypothesis (1.3), this corollary implies that the 
terms in the series in (0.8) will vanish for J^Li ai > deg w (/i). 

Proof of (1.17): We begin by noting that, as in the case with unit weights [25, IV. 20.1]: 

(1.19) If P and Q are weighted homogeneous polynomials in C[x], and deg w (P) — 

-P(x) dx 

deg w (<5) + Yl w i 7^ > then the form u = is exact. 

Q{x) 

Indeed, we find that uj = (deg w (P) — deg w (Q) + Xl^) -1 dcr , where 

P(x) n 

— — — > (— lY~ 1 Wn Xj dxi A ... A dxj A ... A dx n . 

Q(x) 



a = 



The verification of this equality is a straightforward consequence of Euler's formula for 
weighted homogeneous polynomials: 

n dP 
H w i x 3 = deg w (P)P. 



dxj 

j= i J 



As in Lemma (1.12), we write 



R**(t) = J2 ( I x a S m (x)dx ] f 

m> vnv J 
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Since, by Lemma (1.15), B m {x) is a quotient of weighted homogeneous polynomials we 
can apply (1.19) to conclude that 

/ x a B m (x)dx = 

Mk) 

whenever 

deg w (x a ) + deg w (S m ) + w t ^ . 

This inequation is equivalent tom^ s(a). Hence all integrals in (1.12) vanish, except for 
the one with m = s(a). This was precisely the claim of (1.17). o 

The second assertion of (1.17) says that we may write 

(120) Res ( - -J— f ^(f ( - >+1 WB, (a) (x)) 

The numerator is a weighted homogeneous polynomial, by Lemma (1.15). Therefore (1.20) 
is a residue with respect to the (s(a) + 1) power of the initial forms pi(x), . . . ,p n (x). 

We conclude this section by observing that as a direct consequence of (1.17) and the 
Duality Theorem (0.7) we obtain (see [25, IV. 20.1] for the case of unit weights): 

(1.21) Macaulay's Theorem. Letpi(x), . . . ,p n (x) be weighted homogeneous polynomi- 
als whose only common zero is the origin. Then, any weighted homogeneous polynomial 
h(x) satisfying 

n 

deg w (/i) > d w - 

i=i 

is in the ideal generated by pi (x) , . . . , p n (x) . 



2. Grobner Bases for a Term Order 

In this section we specialize to the case of n polynomials gi, ■ ■ ■ ,g n € S with finitely many 
roots in C n , which are a Grobner basis with respect to a term order -<. As in (1.8), we 
choose a positive weight w e N n such that in^(^j) = in w (^j). We may assume that 

(2.1) in w (0i) = x? +1 , i=l,...,n. 

Let r = (n, . . . , r n ) and r + 1 = (n + 1, . . . , r n + 1). Then, with notation as in §1, 
d w = (w, r + 1), and 

P(x) =x[ 1+1 ...x r n " +1 =x r+1 , 



13 



s(a) = (w, a) - d w + ^2wi = (w,&- r) , 
and (1.15) implies the following homogeneity of the coefficients of (1.11): 

(2.2) S m (x) is a w -homogeneous Laurent polynomial of weighted degree — (m + <i w ). 

Consequently, the integrand in (1.17) is a Laurent polynomial. Since 

/ x b dx = for b ^ (-1,...,-1), 
we obtain the following result. 

(2.3) Theorem. The residue Res g (x a ) equals the x^ 1 . . . x~ l -coefficient of x a £? s ( a )(x) . 

Theorem (2.3) is essentially a restatement of a formula for computing global residues due to 
Aizenberg and Tsikh [1, (21.5)], which, in a simpler version pointed out to us by J. Petean, 
says that Res g (x a ) equals the -coefficient in the expression 

e ((-i) ia| x a - (r+i) nfe(x)/x- +i ) 

|a|<(w,a-r) V i=l 

The introduction of the homogenizing parameter t organizes the computation of this Lau- 
rent series and the search for the desired coefficient, as evidenced in Algorithm (3.1) below. 
Keeping track of the homogeneity properties of the coefficients Bj (x) , also allows us to get 
more precise information about the global residues, such as Theorems (2.5) and (2.7). 

(2.4) Remark. Note that (2.2) implies that x a £? m (x) may contain a term of the form 
a x^ 1 . . . x~ x only if 

(w, a) - (m + d w ) = - ^ Wi 

that is, only if 

m = (w, a — r). 

Combined with (1.12), this gives a simpler proof of the first statement in (1.17) in the case 
when <7i, . . . , g n are a Grobner basis with respect to some term order. 

A similar argument combined with the recursive relations (1.16) makes it possible to im- 
prove on the vanishing statement in Corollary (1.18). Let W* denote the polyhedral cone 
in M n which is positively spanned by all lattice points of the form r + 1 — b, where x b runs 
over all monomials appearing in the expansion of (7i(x)(7 2 (x) . . . (7 n (x). 

(2.5) Theorem. The residue Res g (x a ) vanishes if a — r lies outside the cone W* . 




14 



Proof: The condition on b in the definition of W* is equivalent to saying that x b appears 
in one of the coefficients A/(x), j = l,...,d w , in the expansion (1.10) of G(t;x) as a 
polynomial in t. The recursion relations (1.16) imply that B m {x) consists of terms /c Q x a 
where the n-tuples a G Z n are in the translated cone — ((r + 1) + W*). Indeed, if /c a x a is 
a term in B m (x), then (1.16) implies that for some j = 1, . . . ,m, the Laurent polynomial 
Aj(x) S m _j(x) contains a term of the form c a x a+r+1 and therefore, 

a + r + 1 = b + 0, 

where x b is a monomial in Aj{x) and x' 3 is a monomial in S m _j(x). Then 

a + r + 1 = (b - (r+ 1)) + (/? + r + 1) 

and the assertion follows by induction on m. 

Now, Theorem (2.3) implies that Res g (x a ) = unless the coefficient -B s(a )(x) contains 
a term which is a non-zero multiple of x^ ai+r> . . But this is possible only if 

-(a + 1) G -((r + 1) + W*), or equivalently, if a — r G W*. o 



Theorem (2.3) may also be used to study the dependence of Res g (x a ) on the coefficients 
of the polynomials g\, . . . ,g n . We write 

(2.6) ^(x)=^ +1 -f> jX ^, 

i=i 

and let W C R n denote the closed convex cone of all vectors w G R n such that 

(w, aij) < (w, (n + l)ei) ; for all i = 1 . . . , n ; j = 1, . . . , I/*. 

This cone is the polar dual of the cone W* defined above. By assumption, W has non- 
empty interior. Note that W is the cone in the Grobner fan of I corresponding to the 
given term order (see e.g. [14, §3.1]). We have the following result. 

(2.7) Theorem. The residue Res g (x a ) is a polynomial function in the coefficients C{j. 
Its degree in the variable is bounded above by 

(w, a — r) 

(2.8) min 



wew (w, (n + 1) e; - aij) 
and the total degree in the variables c = (cij) is bounded by 

(2.9) min (w, a — r). 



15 



Proof: Given #i(x) as in (2.6), its weighted homogenization with respect to a weight 
w = (wi, . . . , w n ) G Z n fl W , is given by 

x) = - Cij t {w '^ +1) e< - a< ^ x ai ^ 

i=i 

where denotes the i-th unit vector. Set pij = (r* + 1) — a^-. 

According to (2.3), Res g (x a ) equals the ^coefficient of x a £? s ( a )(x), which is 

equal to the coefficient of t< w ' a_r > x^ 1 . . .x~ x in the expansion of 

(2.10) (£ ci,^ w '^>x-^ ) ... £c ni ^)x-^ 

ii,...,i„>0 \J = 1 / \ J = 1 

It is now clear that the residue depends polynomially on the coefficients c^- and that, for a 
given choice of w 6 W, its degree in c^- is bounded above by (w, a— r)/(w, (r^+l) e*— a^-) . 

To prove the bound in (2.9) it suffices to observe that if a monomial c k appears in Res g (x a ), 
then by (2.10) 

J2hj{w,Pij) = (w,a-r). 

i,3 

As (w, p^) > 1 for w G mt(W) integral and all z, j, the claim follows, o 

(2.11) Remarks : i) By the Cauchy-Schwarz inequality, the bound (2.8) is minimized 
by vectors w G W which are as far as possible from a — r and as close as possible to 
(n + l)e t - snj. 

ii) If we are given a system of polynomials /i, . . . , f n which satisfy (2.1), and supposing 
that only the constant coefficients are perturbed, say 



(/i(x) = /i(x)-Ci i = l, ...,n, 

then the bound in (2.9) implies the following: If a monomial c^ 1 . . . cjj™ appears in Res g (x a ), 
then (w, k) < (w, a — r), that is, the weighted degree with respect to w of Res g (x a ) in 
(ci, . . . , c n ) is bounded by (w, a — r) for all integral vectors w G int{W). 



3. Deformation Algorithms for Global Residues 

Let {<7i, . . . ,g n } C S be a Grobner basis as in §2, and let w G N n such that (2.1) holds. 
We have shown that for any polynomial 

h = c a x a G S, 

(w,a)<d 
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the computation of the global residue 



Res 




91 ■ ■ ■ 9n 



x a dx. 



can be reduced to a sum of residues (at the origin) with respect to the family of monomials 
{ . . . , } and some of their powers. Theorem (2.3) gives the following algorithm 

for computing all global residues up to a given weighted degree. Note that Algorithm (3.1) 
respects possible sparsity of the input polynomials. 

(3.1) Algorithm. 

Input: w G W 1 , <?i, . . . , g n G S satisfying (2.1), and deN. 

Output: The global residues Res g (x a ), for all aeN" such that (w, a) < d. 

Step 1: Define the weighted homogenizations (ji(t; x) = t Wi ^ n+1 > gi(t~ Wl xi, . . . , t~ Wn x n ). 
Step 2: If (w, a - r) < 0, then Res g (x a ) = 0. 

Step 3: Set d' = d — (w,r). Compute the Taylor polynomial J2j =0 Bj(x) t J of degree d! 



Laurent polynomial -B( Wja -r) ( x )- R equals Res g (x a ). 

The verification that Q = {<7i, . . . , g n } is a Grobner basis and, if so, the choice of a compat- 
ible weight w G N n may be accomplished in at most 0(n n+2 m 2ri ~ 1 <i( 2n + 1 ) n ) arithmetic 
operations, where m is a bound for the number of monomials in each g^ and d is a bound 
for their degrees. This was shown in [14, §3.2]. 

Naturally, Algorithm (3.1) will be most efficient when it is known a priori that Q is a 
Grobner basis and a "small" compatible weight is given. However, even if the given equa- 
tions Q are not a Grobner basis for any w, then Algorithm (3.1) still serves as a useful 
subroutine. To illustrate this, we describe a general procedure for computing the global 
residue associated to any complete intersection zero dimensional ideal: 

(3.2) Algorithm. 

Input: Polynomials g±, . . . , g n in S whose ideal I is zero-dimensional. 
Output: The global residue Res g (/i), for any specified polynomial h G S. 

Step 1: Choose a "good" term order "-<". 

Step 2: Starting with {<7i, . . . , g n }, run the Buchberger algorithm towards a Grobner 

basis, until the current basis of I contains polynomials fi,...,f n with in^(/j) = 

Step 3: By keeping track of coefficients during Step 2, we obtain an n x n-matrix A = (Aij) 
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of polynomials such that fi = Yl7=i Aij 9j^ or i = 1, ■ ■ • ,n. 
Step 4: Compute the desired residue via the following formula: 

„ / hdx \ ^ fhdet(A)dx 

(3.3) Res = Res ' 



gi---g n J \ h ■ ■ ■ /i 



(3.4) Remarks, (i) In Step 1 we may take "-<" to be optimal in the precise sense of [14, 
§3.3]. Such a choice is possible at almost no extra computational cost if the Grobner basis 
detection procedure of [14, §3.2] had been run beforehand to test applicability of (3.1). 

(ii) The termination and correctness of Step 2 follows from dim(in^(/)) = dim(J) = 0. 

(iii) The correctness of (3.3) is just the Tranformation Law (0.6). In order the evaluate 
the right hand side of (3.3), we may use either Algorithm (3.1), in case many residues are 
desired, or the formula to be presented in (4.2) below, in case only one residue is desired. 

(iv) As shown in [9, Theorem 3.3], it is possible to find polynomials /i, . . . , f n and A^, 
1 < hj < n with degrees bounded by nd 2n + d n + d, where d = max(3, max{deg((7j)}). 

(v) If the polynomials g±, . . . ,g n satisfy (1.3) relative to some weight w, but not necessarily 
(2.1), then the weighted version of the following argument due to Tsikh [25, II. 8. 3] describes 
how to find polynomials fi, ■ ■ ■ , f n in the ideal I = (gi, . . . , g n ) such that 

(3.5) in w (/j) = x? for some p. 

Indeed, Macaulay's Theorem (1.21) implies that, if p > d w — (w± + . . . + w n ), then 

xl G in w (J) = (pi,...,p„). 

This degree bound allows the use of linear algebra (over K) to determine polynomials A^, 
1 < i, j < n, such that = J2i=i A ijPj- Moreover, deg w (A ij ) = pwi - deg w (^). Let 

n n 

fi ■= ^2 9j = $1 + ^2 Aij Qj 

i=l i=l 

Since deg w (Yl A^ qj) < pwi, (3.5) holds, and (3.1) or (4.2) are applicable. Naturally, we 
may use the Buchberger Algorithm to organize the computation of the A^j, which leads 
to a version of Algorithm (3.2) which is applied to the initial forms p±, . . . ,p n rather than 
the entire equations g±, . . . , g n . In summary, there is plenty of room for experimentation ! 

The fact that the coefficients Bj (x) are weighted-homogeneous implies that we can fix the 
value of one of the variables, say x n = 1, and obtain a non- homogeneous version of Algo- 
rithm (3.1): Set x n = 1 everywhere, and, for each a with (w, a) < d, find the coefficient 
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of Xi ' . . .x n _ 1 ™" 1 in -B( w ,a-r)(2 ; i7 • • • >^n-i) !)• We note that other terms of the 
form kx 1 ( ai+1 ) . . . x n l a 1 ™~ 1+1 ' ) may appear in different coefficients Bj(xi, . . . , x n _i, 1), but 
the homogenizing parameter t keeps track of the only one contributing to the residue. 

In the classical case w = (1, . . . , 1), we can apply an argument of Yuzhakov [27] to give 

the following geometric interpretation of (3.1). As in §1, we imbed C n in P n ; the n-form 
x a <ix 

— — — — — may be extended to a meromorphic form $ in P n , which has {t = 0} as a 

polar divisor of order (w, a — r) + 1 if and only if (w, a — r) > 0. The global residue in C n 
may now be expressed as a single residue at a point at oo: if P = (0, . . . , 0, 1), then 

Res g (x a ) = (-l) n Res P $ = 

1 / <9< w > a - r > / i \ \ 

ryReso x? 1 . . . afc 1 • ow w a _ r) ffn — — ( 7-t ) (0) dx\ A ... A dXn-A 

(w,a-r)! V 9t { ' r) \Ui=i9i( x u--^ x n-iA\t)J J 

which, in turn, may be seen to equal the coefficient of x\ ( ^ ai+1 - ) . . . £ n i a 1 "~ 1+1 ' ) in 
-S(w,a-r) ( x ij ■ ■ ■ ,x n -i, 1). With suitable modifications, the same interpretation holds for 
arbitrary weights. 

4. Normal Forms 

In this section we give a formula expressing the coefficients of the (Grobner basis) normal 
form in terms of global residues. In particular, the global residue of a polynomial equals 
the highest coefficient of its normal form (4.2). This leads to a fast algorithm for computing 
residues as well as traces over a zero-dimensional complete intersection. 

Suppose gi,...,g n G K[x] satisfy (2.1) with in w (^) = x^ i+1 . Then V = K[x]/I is an 
Artinian ring of K-dimension (n + 1)(?"2 + 1) • • • (r n + 1). Abbreviating 

J := {i= (ii,...,i n ) eZ™ : < ij < rj , j = 1, . . . , n }, 

the set of monomials {x 1 : i G 1} is a K- vector space basis of V. Every polynomial 
ft 6 K[x] has a unique normal form 

(4.1) NF{h) = J^Ci^x 1 . 

iei 

The scalars Ci(ft) G K are uniquely defined by the property that h = J\fJ-(h) (mod I). 
They are computed using the division algorithm modulo the Grobner basis {gi, . . . , g n }. 

(4.2) Lemma. With the notation as above, every polynomial h G K[x] satisfies 

/ ftrfx \ 

Res — = c ri,...,r„(ft) 

V 9l ■ ■ ■ 9n J 
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Proof: By linearity of the residue operator, Res g (/i) = Xaex CiRes g (x 1 ). However, for 
i G I, i 7^ r, we have (w,i — r) < and, consequently, Res g (x 1 ) = by Corollary (1.18). 
On the other hand, it follows from Theorem (1.17) that 



Res 



x r <ix 



.01 ■■■9n 
This proves Lemma (4.2). o 



Res 



x r dx. 



!"!+ 1 ^n + l 



X 



. . . "JCj 



We now show that all coefficients of the normal form may be computed using residues: 

(4.3) Theorem. Fix an order on the index set X, and let M be the symmetric \T\ x \T\ 
matrix defined by 

x 1 x J dx \ . . 

, i,jeJ. 



My := Res 
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Then M is invertible and for any h G K[x], 

(c i(ft >W = M _1 ( Res (~~~)) J 



ex 



Proof: The Duality Theorem (0.7) implies that the symmetric bilinear form 

hi fi2 diC 



(4.4) 



V xV — > K , (/ii,/i 2 ) ^ Res 
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is non-degenerate. The symmetric matrix M represents (4.4) relative to the basis { x 1 
i El}. Therefore M is non-singular. The second claim follows from the fact that 



Res 



x 1 h <ix 
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= J2cj(h)Res(^^-) = £M„ C j(fc).o 
jex \gi...g n J . £J 



(4.5) Remark. As observed in [3], if we choose the lexicographical order in J, then the 
matrix M has the triangular form 



M 



( 



Vi 



1\ 



Consequently, detM = ±1 and it is easy to compute the inverse of M. 
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We have seen in §0 that the Duality Law may also be interpreted as saying that the global 
residue operator Res g is a generator of V* as a F-module. There is another element in 
this F-module which is of special interest in computational algebra, namely, the morphism 
tr G V* which assigns to a polynomial h the trace of the endomorphism of V given by 
multiplication by h. The trace can be computed by normal form reduction as follows: 

(4.6) tr(h) = ^Ci^-x 1 ). 

iei 

On the other hand, it is known (see e.g. [23, Satz (4.2)]) that the trace may be expressed 
in terms of the global residue: 

(4.7) tr(h) = Yl VMHP) = Res 

pez( s ) 

This expression, together with (4.2), gives the following formula for computing the trace: 

(4.8) Algorithm. Compute the trace of an element h G V by tr(/i) = c ri ___ rn (h ■ J g ). 

Thus, to find the trace of an element of V over K, it suffices to run a single normal form 
reduction. We found Algorithm (4.8) to be quite efficient in practice. Additional speed 
can be gained by simple tricks, such as replacing J g by AfJ r (J s ) in Algorithm (4.8), and 
by storing previously computed normal forms of monomials. 

From Theorem 2.7 we can derive bounds on the degree of the trace for gi, . . . , g n as in (2.6) 
which satisfy (2.1). Note that the corresponding cone W in the Grobner fan has nonempty 
interior. For each value of the parameters c^-, let Z c denote the zero set of g\, . . . ,g n . 

(4.9) Theorem. For any h G K\x\, the parametric trace 

tr(/i)(c) = ^^(pJMp) 

pez c 

is a polynomial function of c = (q 7 -) with degree bounded above by min deg w (/i). 

w6int(W)nZ™ 

Proof : Let J g (x) be the Jacobian of g\, . . . ,g n with respect to the variables x±, . . . , x n . 
Given any polynomial h and w G int(W) fl Z n , we know by the second bound in The- 
orem 2.7 that the degree of tr(/i)(c) = Res g (/i • J g ) in the c variables is bounded by 
deg w (/i • J g ) — (w,r). As deg w (J g ) = (w,r), the claim follows, o 

5. Real Roots, Degree, and Symmetric Polynomials 



h J g <ix 
9i ■ ■ ■ 9n 
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We present three applications of the computation of global residues: counting real roots 
using the trace form (following Pedersen-Roy-Szpirglas [22], Becker- Wormann [4]), com- 
puting the degree of a polynomial map (following Eisenbud-Levine [12]), and evaluating 
elementary symmetric polynomials in a multivariate setting (following classical work of 
Junker [16]). 

We assume as above that g±, . . . , g n E K[x] satisfy (2.1) with in w (<7j) = x^ i+ , and again 
let X := { i = (ii, . . . , i n ) E Z n : < ij < Tj , j = 1, . . . , n }. Suppose that K is a 
subfield of the real numbers M and let T be the symmetric |X| x |X|-matrix T defined over 
Kby 

(5.1) Tij := tr^xi), ijGl. 

The following result is due to Becker- Wormann and Pedersen-Roy-Szpirglas. 

(5.2) Theorem. ([4], [22]) The rank ofT equals the number of distinct complex roots 
in Z(g). The signature ofT equals the number of distinct real roots in Z(g) D M. n . 

Recall that the signature of T equals the number of positive eigenvalues minus the number 
of negative eigenvalues (all eigenvalues of a real symmetric matrix are real). As is pointed 
out in [22, Prop. 2.8], the signature can be read off directly from (the number of sign 
variations in) the characteristic polynomial of T. A straightforward generalization of (5.2) 
states that, for any h E V, the signature of the matrix T h = (tr(x'x j ■ h)) i . eJ equals 
the number of distinct real roots with h > minus the number of distinct real roots with 
h < 0. Algorithms (4.8) and (3.1) provide subroutines for computing T and hence for 
counting real zeros of zero-dimensional complete intersections. 

Viewing now (<7i, . . . , g n ) as a proper map g : R n — > M n , its degree is defined as 

deg(g) := Yl de S P (g) 

where q is a regular value of g and deg p (g) is ±1 depending on whether J s (p) is positive 
or negative. The degree is a topological invariant of g. 

Let M be the non-singular, symmetric matrix M defined, as in (4.3), by 

^ / x 1 x^ dx \ . . _ 

My := Res , i,j El. 

\9l---9nJ 

The following result is essentially contained in [12]; although the results there are local, 
the passage to the global situation may be done as in [22]. 

(5.3) Theorem. The degree of g equals the signature of M. 
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We may apply Algorithm (3.1) or Lemma (4.2) to compute the matrix M and, conse- 
quently, the degree of g. 

For our third application we need to review the concept of symmetric polynomials in a 
multivariate setting. This theory is classical (see Junker [16], who refers to even earlier 
work of MacMahon and Schlafli). It reappeared in the recent computer algebra literature 
in [21]. Let A = (otij) be an N x n-matrix of indeterminates over K. The symmetric group 
Sn acts on the polynomial ring KfcKy] by permuting rows of A. We are interested in the 
invariant subring K.[aij] SN , whose elements are called symmetric polynomials. It is known 
that K.[ctij] SN is generated by symmetric polynomials of total degree at most N, but, in 
contrast to the familiar n = 1 case, this K-algebra K[ajj] Sjv is not free for n > 2. An 
important set of generators are the elementary symmetric polynomials e$(A), which are 
defined as the coefficients of the following auxiliary polynomial in m, 112, ■ ■ - u n : 

N 

(5-4) Y[(l + a il u 1 + a i 2U2-\ hft m tl„) = ^ e ji,...,j n ( A ) ' u l u 2 • ' • u n 

i=l ji+-+j n <N 

Another set of generators is given by the power sums : 

N 

(5.5) h s := J2< a %---<' for i = (hJ2,...Jn) e N», ji + ... + j n <N 

1=1 

Algorithms and formulas for writing the e- 3 in terms of the h- 3 and conversely are studied 
in detail by Junker [16]. One of his methods will be presented in (5.8)-(5.9) below. 

Returning to our zero-dimensional complete intersection, let = dimK(V^) be the cardi- 
nality of the multiset Z(g) C C n (counting multiplicities). We fix any bijection between 
the rows of A = (a^) and Z(g). This defines a natural K-algebra homomorphism 

(5.6) <j> : K[a>ij] SN ^ K, 

where the indeterminate gets mapped to the j-th coordinate of the i-th point in Z(g). 
Our objective is to evaluate the map using only operations in K. In particular, we are 
interested in the problem of evaluating the elementary symmetric polynomials ej under <p. 

The punch line of our discussion is that it is easy to evaluate the power sums via the trace: 

(5.7) = tr(xJ). 

Thus to compute (5.7) we use Algorithm (4.8). We then proceed using the following 
method due to Junker and MacMahon. Consider the image of (5.4) under 0, 

(5.8) R(u) = H (l+p lUl + ...+ PnUn )*««(p) = ^0( ej ). U J 

pez( s ) j 
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The polynomial R(u) is the Chow form of the zero-dimensional scheme defined by /. In 
computer algebra it is also known as the U-resultant. Following[16, pp. 233, Eq. (4)], the 
formal logarithm of (5.8) equals 

(5.9) io g (i*(u)) = £ tn^l . y: (f) m) 

d=l a \i\=d vjy 

Here |j| = ji + ■ ■ ■ + j n an d (j) = j^t'-.j Using (5.7) and (4.8), we can compute the 
formal power series (5.9) up to any desired degree 61 '. We then formally exponentiate this 
truncated series (using operations only in K) to get the Chow form (5.8) up to the same 
degree 6'. In order to determine (5.8) completely, which means to evaluate all elementary 
symmetric polynomials, it suffices to expand (5.9) up to degree 61 = N = diniK(U). 



6. An Example 

In this section we apply our results and algorithms to the specific trivariate system: 

(6.1) g\ = x\ + x\ + x\ - 1, g 2 = x\ + x% + x 3 - 1, g 3 = x\ + x\ + x\ - 1. 

This example is taken from [14, Example 3.1.2], where it served to illustrate the problem of 
Grobner Basis Detection. Indeed, the polynomials <7i,<72,<73 are a Grobner basis, namely, 
for the weight vector w = (3,4, 7). With respect to these weights, the initial monomials 
are the pure powers underlined above. We see that, counting possible multiplicities, the 
set Z(g) consists of 30 points in C 3 . Our basic problem is to evaluate the global residue 

(6.2) Res g (x a ) = Res ' Xl * 2 ^ * 



,#l(x)#2(x)#3(x) 

for any non-negative integer vector a = (ai, 0,2, 0,3). 

It is interesting to compare the relative efficiency of Algorithm (3.1) and the Grobner basis 
reduction method deduced from Lemma (4.2). In step 1 of Algorithm (3.1) we compute 
the weighted homogenizations 

g x = x \ + t 3 xl + t 7 xl -t 15 , g 2 = x\ +tx 3 + t 2 xl -t 8 , g 3 = x\ + tx^ + t 3 x 6 1 -t 21 . 
We then consider the expression 

1 



(6.3) 



(?i(£;x) • g 2 (t;x) ■ g 3 (t;x) 
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as a rational function in t, and we compute its Taylor expansion Y^=o Bj ( x ) ^ U P to some 
degree d which exceeds 3(ai — 4) + 4(a 2 — 1) + 7(a 3 — 2). Here the coefficients Bj(x) are 
w-homogeneous Laurent polynomials in xi,X2,X3; for instance, 

_ 1 1 x-2 x\ ./•;{ x\ 1 

2\-X-J Jo 4 3" 4" 3 5 5 10 4 15 2 5 9 5 6 ' 

^1 *^2 X^XqX^ X-^X3 X-^ ^3 X-^ ^2 X-^X3 tC^^2^3 

Now set j = 3(ai— 4)+4(a2 — l) + 7(a3 — 2). The desired residue (6.2) equals the coefficient 
of Xi ai ~ 1 x 2 ~ a2 ~ 1 X3 0-3-1 in the Laurent polynomial Bj e Z[xi, xj -1 , x 2 , x^ 1 , X3, X3 1 ]. 

This Taylor expansion is a fairly space consuming process since the polynomials Bj(x) 
grow quite large. This is witnessed by the following table, which shows the number of 
terms of Bj(x) for some values of j between 2 and 40: 

j : 2 5 10 15 20 25 30 35 40 
# : 7 41 216 569 1102 1803 2682 3744 4964 

On the other hand, the normal form method of Lemma (4.2) is quite efficient for evaluating 
individual residues. Let I denote the ideal in Q[xi, X2, X3] generated by (6.1). The quotient 
ring V = Q[xi,X2,X3]/I is a 30-dimensional Q-vector space. Every element h G V is 
uniquely represented by its normal form NJ-{h) modulo the reduction relations: 

(6.4) x\ — ► -x\ - x\ + 1, x\ — > -xj-x 3 + l, x\ — > -xf-xl + l. 

By Lemma (4.2), the residue (6.2) is equal to the coefficient of x\x2x\ in MT{K). 
For instance, for the Jacobian 

J(x) = det (— -) = 18x^X2 — 2Ax\x 2 x% — 2hx\x% + 30xfx2X^ + 2OX1X2X3 — ld>x\x\x\ 

it takes 10 reductions modulo (6.4) to reach the normal form 

AfJ-'(J) = 30x^X2^3 — 2bxfx% — 152x\x 2 + 146x^X3 — 251xfx2X3 + 83xfx| + 16xf 
+ 229x?x 2 + 8x^x3 - 196x?x 2 x 3 + 226xjxl - 114x^x1 - 73x? + 240xjx 2 
+ 34xfx3 + 254xix 2 xs — 62x1X3 + 69^2X3 — 260x^ — 140xiX2 — 78x1X3 
+ 108x 2 x 3 - 49xg + 140xi - 177x2 - 128x 3 + 177. 

Indeed, we see that the coefficient of x\x 2 x\ equals Res g (J) = tr(l) = dim (V) = 30. 
Here is a slightly more serious example: it takes 62 reductions modulo (6.4), running less 
than two minutes in MAPLE on a Sparc 2, in order to find the global residue 

Resg(x} 5 x^ 5 ) = -258,756,707,658,424,020,014,953,731,203. 
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We made the observation that the efficiency of the two methods is comparable when 
computing all residues of the form Res g (x a ) with (w, a) < d for some fixed d. This is the 
case, for example, in the computation of the matrix M defined in §4. This is a symmetric, 
30 x 30-matrix whose computation using Algorithm (3.1) requires the knowledge of -Bj(x) 
for j < 30. Using MAPLE on a Sparc 2 these may be obtained in 321 seconds. It takes 
an additional 247 seconds to read off the desired 465 coefficients. On the other hand, it 
takes 324 seconds to build up the matrix M using Lemma (4.2). The signature of M is 
zero, and hence so is the degree of the map g : R 3 — > R 3 by Theorem (5.3). 

For further combinatorial analysis we may wish to compute the two polyhedral cones in 
Section 2. We first obtain the 3-dimensional quadrangular cone 

W = { (wi, W2, ws) G R 3 : 5tui > 2ty 3 , 2ty 2 > ^3, ^3 > 2tui, 3-103 > 5w 2 } 
= pos{(4,5,10), (1,1,2), (5,6,10), (2,3,5)}. 

The interior of W consists of all weight vectors which select the underlined monomials in 
(6.1) to be initial. The cone polar to W equals 

W* = pos {(5, 0,-2), (0,2,-1), (-2,0,1), (0,-5,3)} 

= { (ai, <22, 03) G -R 3 : 4ai + 5a2 + 10a3 > 0, ai + a 2 + 2a3 > 0, 

5ai + 6a 2 + 10a 3 > 0, 2ai + 3a 2 + 5a 3 > 0}. 

By Theorem (2.5), the residue (6.2) vanishes whenever 

(ai, a2, as) G" (4, 1, 2) + W*, or equivalently, 

4ai+5a2 + 10a3 < 41 or ai+a2+2a3 < 9 or 5ai+6a2 + 10a3 < 46 or 2ai+3a2+5a3 < 21. 
For instance, (6, 1, 1) satisfies the first inequality and therefore Res g (a;f X2X3) = 0. 

In Section 4 we have shown that the trace tr(/i) of an element h in V = Q[x]/7 can be 
computed easily as the coefficient of x\x2x\ in MT{h- J). Using this technique, let us now 
analyze the zero set Z(g) with respect to multiple roots, real roots, etc... We compute the 
symmetric, integer 30 x 30-matrix representing the trace form T as in (4.8). The largest 
entry in T appears in the lower right corner: 

tr{x\x 2 xl ■ x\x 2 x\) = Resg^fa;^! • J(x)) = 16,049,138,278. 

The rank of the matrix T equals 20. By Theorem (5.2), this is the number of distinct roots 
of g. The characteristic polynomial of T has 13 positive real roots and 7 negative real 
roots. Therefore the signature of T equals 6, and this is the number of distinct real roots 
of g. It turns out that there are four rational roots, and they account for all multiplicities: 
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the root (1,0,0) has multiplicity 3, the root (0,1,0) has multiplicity 4, the root (0,0,1) 
has multiplicity 6, while the root (—1, 1, —1) is simple. The remaining 16 roots, two real 
and 14 imaginary, are all simple and they are conjugates over Q. 

We finally come to the problem of computing the Chow form 

R{u u u2,uz) = Y\ (1 + aiwi + a 2 «2 + a3« 3 ) Atg(a) - 

Note that each of the three non-simple roots appears with its multiplicity in this product. 
The ( 3 3 3 ) = 5,456 rational coefficients of R(u\, 112, 113) are the values of the elementary 
symmetric polynomials at the roots of g = (g±, g2, (73). Following (5.9), (5.7) and using 
Algorithm (4.8), we compute the following formal power series up to a chosen degree: 

log(i?(-Ui, U2, u 3 )) = tr(xi)-ui + tr(x 2 )u 2 + tr(x 3 )u 3 - ^ • (tr(a^)^ + 2tr(x 1 x 2 )u 1 U2 
+ 2tr(x 1 x 3 )u 1 us + tr(x\)u\ + 2tr(x 2 x 3 )u 2 us + tr(x§)u|) + ^ • (tr(xf)-u 3 + . . . 

O 

35 2 485 2 o o 

= 5ii2 — 5ii3 + 37uiit,2 — 121uiii3 — — u 2 + IO6W2W3 — u s + 17u 1 — 7\u x u 2 

+ mu-^us — YJ2u\u 2 + 536witt2^3 _ 686wiw 3 H — — u 2 — 667u 2 us + 1084u2U 3 + . . . 
By formally exponentiating this series, we obtain the Chow form 

R(u±, k.2, U3) = 1 + 5-U2 — 5tt3 + 37uiU2 — 121-U1-U3 — 5u 2 + 81^2^3 — 230tt3 + 17-uf 
— 7Au\u2 + 177u1us + 13uiu 2 — 254uiU2U3 — 81uiu 3 — 5u 2 — 112u 2 us — 596^2^3 + . . . 

and hence all elementary symmetric polynomials. For instance, we see that X^ a i/^273 = 
—254, where the sum is taken over all triples of roots (a±, «2, 0^3), (Pi, $2, Z^), (71, 72, 73) £ 
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Note added in proof: During the MEGA 94 meeting we became aware of the paper: [M. 
Kreuzer and E. Kunz: Traces in strict Frobenius algebras and strict complete intersections. 
J. reine angew. Math. 381 (1987), 181-204]. Our assumption (1.3) is equivalent, by their 
Proposition (4.2), to the statement that the iiT-algebra V = K\x\/I is a strict complete 
intersection. Consequently, Theorem (1.17), in the case s(a) < (in particular the Euler- 
Jacobi Theorem (Corollary (1.18)) is contained in their Corollary (4.6) and Theorem (4.8) 
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